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An  Outline  for  the  RTNEPH  Infrared 
Atmospheric  Water  Vapor  Attenuation  Study 


1.  BACKGROUND 

Real-Time  Nephanalysis  (RTNEPH)  is  an  automated  cloud  analysis  model  that  is 
in  operational  use  at  the  Air  Force  Global  Weather  Central  (AFGWC).  In  the 
RTNEPH,  polar-orbiting-satellite  imagery  is  analyzed  in  conjunction  with 
conventional  meteorological  cloud  observations  to  produce  a  global  analysis  of  cloud 
attributes  such  as  extent,  height,  bases,  and  type.  The  RTNEPH  and  its  predecessor, 
the  Three-Dimensional  Nephanalysis  (3DNEPH),  have  been  generating  real-time 
global  cloud  analyses  since  1970  (see  Fye,  1978;  Kiess  and  Cox,  1988).  RTNEPH 
cloud  analyses  are  used  primarily  to  initialize  cloud  forecast  models  in  support  of  Air 
Force  missions;  they  are  also  archived  for  use  by  the  research  community  as  one  of 
the  few  sources  of  long-term  global  cloud  climatologies  available  today. 

RTNEPH  makes  extensive  use  of  satellite  data  measured  in  the  10-12p.m  infrared 
(IR)  window  for  the  thermal  mapping  of  clouds  and  the  earth’s  surface  both  day  and 
night.  Although  at  first  approximation  this  part  of  the  spectrum  is  generally 
considered  free  of  water  vapor  absorption  effects,  attenuation  is  significant  for  high 
viewing  angles  and  tropical  atmospheres.  Thermal  radiation  emitted  by  the  ground  is 
attenuated,  primarily  by  water  vapor,  as  it  travels  up  through  a  cloud-free 
atmosphere.  The  net  effect  is  that  the  radiance  observed  by  the  satellite  is  less  than 
the  radiance  emitted  by  the  underlying  surface,  causing  that  surface  to  appear  colder 
than  it  truly  is. 

The  RTNEPH  uses  a  threshold  technique  whereby  a  satellite-measured  IR 
brightness  temperature  is  compared  to  the  actual  skin  temperature  of  the  surface 
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being  observed.  If  the  brightness  temperature  is  equal  to  the  true  ground 
temperature,  the  scene  is  considered  cloud  free;  if  it  is  lower,  the  scene  is  determined 
to  contain  cloud.  The  underlying  assumptions  are  that:  1)  temperature  generally 
decreases  with  height  throughout  the  troposphere,  and  2)  clouds  are  always  physically 
colder  than  the  ground;  therefore  their  temperatures  are  always  lower  than  the 
ground  temperature.  (Obviously  the  presence  of  low-level  temperature  inversions 
invalidates  this  assumption.  How  to  handle  inversions  is  an  issue  (aside  from 
attenuation)  that  needs  to  be  addressed  more  carefully  in  the  RTNEPH.) 

Water  vapor  attenuates  the  emitted  ground  radiance,  reducing  its  magnitude  and 
with  it  the  apparent  temperature  of  the  surface  being  viewed.  Without  an 
attenuation  correction  being  applied  to  the  satellite  brightness  temperature 
measurement,  it  will  almost  always  be  lower  than  the  skin  temperature  of  the  surface 
being  viewed.  If  a  raw  satellite-measured  brightness  temperature  is  compared 
directly  with  the  actual  surface  skin  temperature,  the  RTNEPH  will  always  flag  a 
scene  as  containing  clouds.  Thus  it  is  easy  to  see  the  importance  of  properly 
accounting  for  attenuation  losses  in  IR  brightness  temperature  measurements  before 
they  are  processed  by  cloud  detection  algorithms. 

In  recent  years  the  Geophysics  Laboratory  has  been  actively  involved  in 
proposing  and  developing  improvements  for  the  RTNEPH  (Bunting  et  al.,  1983; 
d’Entremont  et  al.,  1982;  d’Entremont  et  al.,  1989).  It  is  the  intent  of  this  report  to 
outline  in  detail  a  new  study  designed  to  enhance  the  cloud  detection  capability  of  the 
RTNEPH  by  more  effectively  handling  the  effects  that  atmospheric  water  vapor 
attenuation  have  on  satellite  brightness  temperature  measurements  at  IR 
wavelengths;  upon  completion,  results  of  this  study  will  be  published  in  a  subsequent 
paper.  First,  atmospheric  water  vapor  attenuation  theory  will  be  presented  in 
Sections  2  and  3,  followed  by  a  discussion  in  Sections  4-7  on  how  to  apply  this  theory 
to  the  RTNEPH  in  as  accurate  and  sensible  a  manner  as  possible.  Section  8  outlines 
the  plans  for  developing  and  incorporating  a  water  vapor  attenuation  model  into  the 
RTNEPH,  taking  into  consideration  both  the  type  and  accuracy  of  the  atmospheric 
data  that  the  RTNEPH  has  (or  soon  will  have)  available  for  use. 
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2.  NOTATION 


Ix  Monochromatic  Upwelling  Thermal  Radiance, 

W  m~2  p  m-1  ster-1  [  M  L-1  T-3  J 

BX(T)  Monochromatic  Planck  Blackbody  Radiance, 

W  m~2  p,m-1  ster-1  [  M  L-1  T-3  ] 

\  Wavelength,  p,m  [  L  ] 

T  Temperature,  K  [  0  ] 

rx(p)  Atmospheric  Transmittance,  Dimensionless 

ex  Monochromatic  Emissivity,  Dimensionless 

rx  Monochromatic  Reflectivity,  Dimensionless 

tx  Monochromatic  Transmissivity,  Dimensionless 

Ij  Average  Spectral  Upwelling  Thermal  Radiance, 

W  m~2  (im'1  ster-1  [  M  L-1  T-3  ] 

4)j(\)  Sensor  Response  Function  for  Satellite  Sensor  j, 

Dimensionless 


3.  RADIATIVE  TRANSFER  THEORY 

The  monochromatic  upwelling  thermal  radiance  at  the  top  of  a  non-scattering, 
plane-parallel,  cloud  free  atmosphere  that  is  in  local  thermodynamic  equilibrium  and 
whose  source  function  is  the  Planck  function  may  be  written  in  the  form 

0  <?T 

h  =  ^,,r,Bx(T,t«Kiltc  +  J  BJT(p')]— 7-dp' 

P.fc 

Pttc  Qyl 

+  r\,*fcT\,*fc  /  Bx[T(p')]— 7-dp',  (1) 

o  dp 

where  eXgfc  is  the  emissivity  of  the  earth’s  surface;  BX(T)  is  the  Planck  radiance 
emitted  by  a  blackbody  of  temperature  T  at  wavelength  TXigfc  is  the  atmospheric 
transmittance  at  wavelength  X.  between  the  top  of  the  atmosphere  and  the  surface; 
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Tx  =  Tx(p)  is  the  atmospheric  transmittance  at  wavelength  X  for  the  vertical  path 
between  the  top  of  the  atmosphere  (at  p=0)  and  any  pressure  level  p  in  the 
atmosphere;  rX3fc  is  the  surface  reflectivity  (rXsfc  =  1— €Xsfc);  and  T'x(p)  is  the 
transmittance  between  the  surface  at  p=psfc  and  any  pressure  level  in  the 
atmosphere  p  (note  that  Tx(psfc)  =  t'x(0)  =  tx  sfc  ).  The  first  term  on  the  right  side 
of  Eq.  (1)  denotes  the  ground  contribution  to  Ix  and  is  the  dominant  term  in  long- 
wavelength  IR  (LWIR)  thermal  window  regions.  The  second  term  denotes  the 
upwelling  atmospheric  emission,  and  the  third  term  denotes  the  atmospheric  emission 
radiated  in  the  downward  direction  and  reflected  back  up  by  the  surface. 

The  Planck  blackbody  emission  BX(T)  for  a  blackbody  of  temperature  T, 
wavelength  A,  is  given  by 


Ba(T) 


2hc2A-5 

ehc/KAT  _  j 


where  A  is  in  meters,  T  is  in  Kelvin,  the  Planck  constant  h  =  6.6262  x  10-34  J  sec, 
the  speed  of  light  c  =  2.998  X  108  m  sec-1,  and  Boltzmann’s  constant  K  =  1.3806  x 
10-23  J  K-1.  For  convenience  it  can  be  rewritten  in  the  form 


Bx(T) 


ctX-5 

e^/xT  -  1  ’ 


(2) 


where  ct  =  1.19107  X  108  W  p,m4  m-2  ster-1,  c2  =  1.43883  x  104  p.m  K,  T  is  in  K, 
and  where  X  is  now  expressed  in  p,m  (instead  of  m).  Using  Eq.  (2),  the  Planck 
radiance  units  are  W  m-2  p.m-1  ster-1. 

Spaceborne  satellite  sensors  are  not  designed  to  measure  monochromatic 
radiances  but  radiant  energies  over  some  wavelength  range  Xj^XsX2.  The  observed 
radiance  Ij  sensed  by  a  downward-pointing  radiometer  in  spectral  band  “j”  is  given 
by 


Xj2 

/  <t>j(*)dX 

^ - ■  (3) 

J  <bj(x)dx 

u. 

where  <j>j(X)  is  the  spectral  response  function  for  sensor  j,  with  the  j  subscript  in  Ij 
denoting  a  spectrally  weighted  average  radiance.  The  response  function  c}>j(M  is  a 
measure  of  sensor  j ’s  instrument  response  to  radiation  at  wavelength  X;  if  cf>j(X)  =  1, 
the  sensor  detects  100  percent  of  the  energy  radiated  at  wavelength  X,  whereas  if 
<t>j(X)  =  0.25,  then  the  sensor  detects  only  25  percent  of  energy  radiated  at 
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wavelength  Within  some  wavelength  range  the  response  function 

cf)j(X)  for  sensor  j  is  a  rapidly  varying  function  of  wavelength;  outside  this  range  it  is 
zero.  Response  functions  are  specific  to  each  sensor  and  indicate  instrument 
sensitivity  characteristics.  Response  functions  for  the  NOAA  Advanced  Very  High 
Resolution  Radiometer  (AVHRR)  thermal  channels  4  and  5  and  the  DMSS 
Operational  Linescan-Thermal  (OLS-T)  channel  are  shown  in  Figure  1. 


4.  APPLICATION  TO  THE  RTNEPH 

Equations  (1)  and  (3)  are  useful  for  computing  water  vapor  attenuation  effects  in 
LWIR  satellite  radiance  measurements.  To  adapt  Eq.  (1)  for  operational  RTNEPH 
use,  the  following  approximations  are  made.  First,  LWIR  emissivities  ex  for  the 
ground  and  most  opaque  clouds  are  assumed  to  be  one.  With  €XjSfc  =  1,  then 
rx,sfc  =  i-ex.sfc  =  0,  and  Eq.  (1)  reduces  to 

°  0T 

R*  =  Bx(T,ft)Tvft  +  /  Bt|T(p')]-rr<ip'.  (4a) 

p.„  3P 

Equation  (4a)  gives  the  form  for  the  monochromatic  radiance  emitted  by  a  blackbody 
surface  (in  this  case,  the  ground)  of  temperature  Tgfc.  Substituting  R*  into  Eq.  (3)  for 
lK  yields 


Xj2 

f  d>j(k)d\ 


f  4>j(Mdx 

xi> 

where  Rj  is  the  satellite-observed  radiance  emitted  by  a  blackbody  surface  and 
travelling  through  a  cloud  free  atmosphere. 

4.1  Cloud/No-Cloud  Decisions 

The  RTNEPH  has  access  to  surface  and  upper  air  data  that  allow  for  the 
calculation  of  the  expected  satellite  radiance  as  given  by  Eq.  (4b).  Surface  skin 
temperatures  Tgfc  are  available  to  the  RTNEPH  on  an  eighth-mesh  (47  km) 
horizontal  grid  resolution.  Thus,  using  Eq.  (2),  the  Planck  radiance  term  Bx(Tgfc)  on 
the  right  side  of  Eq.  (4a)  can  be  computed  on  an  eighth-mesh  basis.  The  temperature 
profile  T(p)  is  available  on  a  whole-mesh  (376  km)  resolution,  allowing  the 
atmospheric  radiances  Bx(T(p))  to  be  computed  on  a  whole-mesh  basis. 
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SENSOR  RESPONSE  (PERCENT) 


figure  1.  Response  Functions  for  the  NOAA  AVHRR  10.7  and  11.8p,m 
LWIR  Channels  4  and  5,  and  the  DMSS  OLS-T  Thermal  Channel 
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The  transmittance  profile  Tx(p)  denotes  the  amount  of  atmo«pheric  attenuation 
at  wavelength  X  between  the  top  of  the  atmosphere  and  the  level  at  pressure  p. 
Transmittances  can  also  be  computed  at  whole-mesh  resolution  since,  in  the  LWIR, 
Tx  is  (essentially)  a  function  of  temperature  and  water  vapor  profiles.  The  details  of 
how  the  RTNEPH  can  compute  transmittances  in  real-time  are  left  for  the  next 
section;  assume  for  now  that  the  Tx(p)  are  provided  at  the  whole-mesh  resolution  of 
the  upper  air  profiles. 

With  surface  skin  temperature  and  the  corresponding  upper  air  temperature  and 
transmittance  profiles,  the  monochromatic  radiance  Rx,  and  in  turn,  the  spectrally 
averaged  radiance  Rj  can  be  computed  for  each  eighth-mesh  grid  point  of  the 
RTNEPH.  Assuming  that  the  atmospheric  temperatures  T(p)  decrease  monotonically 
with  pressure,  there  is  no  way  for  a  satellite-measured  radiance  Isat  from  a  cloudy 
field  of  view  to  be  higher  than  the  clear-column  Rj.  Thus  Rj  can  be  compared 
directly  to  Isat:  if  I9at  <  Rj  ,  the  scene  contains  clouds;  if  Isat  =  Rj  ,  the  scene  is  cloud 
free. 

One  problem  with  using  this  approach  is  that  the  RTNEPH  uses  satellite- 
measured  brightness  temperatures  Tsat  from  the  Satellite  Global  Data  Base  (SGDB) 
to  perform  the  cloud/no-cloud  decision,  not  radiances  Isat.  Thus,  the  computed 
radiances  Rj  must  first  be  converted  to  an  equivalent  brightness  temperature  Tj  using 
a  brightness-temperature-to-radiance  lookup  table  that  is  constructed  in  the  following 
way: 


Xj2 

/  BJT)  <bj(A)dX 

Bj(T)  =  ^7 - ’  <5a) 

where  BX(T)  is  given  by  Eq.  (2).  Equation  (5a)  relates  temperature  T  to  radiance  Bj 
for  satellite  sensor  j.  By  constructing  a  table  of  values  Bj(T)  for  a  range  of 
terrestrially  expected  temperatures  T  (say  190  -  330  K),  it  is  possible  to  map  a 
computed  (or  measured)  radiance  to  a  unique  brightness  temperature  Tj  by 
substituting  Rj  for  Bj  in  Eq.  (5a).  Let  this  relation  be  denoted 

Tj  =  T(R,),  (5b) 

where  Tj  is  the  equivalent  brightness  temperature  corresponding  to  radiance 
measurement  Rj,  as  defined  by  Eq.  (5a). 

Consider  a  scene  in  which  a  satellite  measures  a  particular  brightness 
temperature  Tsat  .  A  decision  can  be  made  as  follows  concerning  whether  the  scene 
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contains  cloud.  First,  determine  the  clear-column  brightness  temperature  Tj  from 
Eq.  (5b)  using  the  theoretically  computed  clear-column  radiance  Rj  from  Eq.  (4). 
Then  use  Tj  to  compute  a  brightness  temperature  difference 

AT,  =  Tsfc  -  Tj  ,  (6) 

where  Tsfc  is  the  surface  skin  temperature  from  the  RTNEPH  database.  When  ATj  is 
added  to  the  observed  temperature  Tsat,  the  result  is  the  true  physical  temperature 
Tsfc  of  the  ground  being  viewed.  A  direct  comparison  can  be  made  between  the 
surface  skin  temperature  Tsfc  and  T*,  where 

T*  =  T.„  +  AT;.  (7) 

If  the  scene  is  cloud  free,  then  T*=Tgfc.  If  clouds  are  present,  then  T*<Tsfc  since  the 
clouds  are  colder;  they  emit  less  radiation,  thereby  showing  a  lower  temperature  than 
the  warmer  ground  below. 

Note  that  the  effects  of  LWIR  atmospheric  attenuation  are  brought  into  the 
cloud/no-cloud  decision  via  the  transmittances  tx  in  the  monochromatic  upwelling 
radiance  Eq.  (4a).  Thus,  assuming  that  the  skin  temperatures  Tsfc  are  reliable,  the 
critical  issue  here  is  to  compute  the  Tx(p)  profiles  as  accurately  as  possible.  This  issue 
is  addressed  in  Section  5. 

4.2  Improved  Cloud  Top  Height  Calculations 

If  a  blackbody  is  placed  at  some  pressure  level  p*  within  the  atmosphere,  then  Rx 
[from  Eq.  (4a)  ]  can  be  written  as  a  function  of  p*  in  the  form 

Rx(p*)  =  Bx(T(p*)1tx(p*)  +  /  Bx(T(p')|-^7dp',  (8a) 

p*  dp 

where  Rx(p*)  is  the  monochromatic  upwelling  radiance  observed  at  the  top  of  a 
cloud-free  atmosphere  within  which  lies  a  blackbody  at  pressure  level  p  =  p*  [which 
can  be  converted  to  an  altitude  z  =  z*  via  the  pressure-height  profile  p(z)  ]. 
Substituting  Rx(p*)  for  Rx  in  (4b)  yields  the  expected  upwelling  radiance  for  spectral 
band  j : 


Xj2 

j  R\(p*)  <J>j(k)d\ 

Rj(p*)  =  ^ - ,  (8b) 

f  <J>j(X)d\ 

N' 
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where  Rj(p*)  is  the  radiance  emitted  through  a  cloud  free  atmosphere  by  a  blackbody 
at  pressure  level  p*,  and  measured  by  satellite  sensor  j.  A  profile  of  expected 
brightness  temperatures  Tj(p*)  can  be  computed  for  different  p*’s  using  Eq.  (8), 
where 


Tj(p*)  =  T[Rj(p*)] 


(8c) 


from  Eq.  (5b). 

Now  consider  the  blackbody  at  pressure  level  p*  to  be  a  cloud.  Assuming  (as  the 
RTNEPH  does)  that  the  cloud  completely  fills  the  pixel  field  of  view,  a  profile  of 
ATj’s  can  be  computed: 


ATj(p)  =  T(p)  -  Tj(p)  ,  (9) 

where  T(p)  is  the  atmospheric  temperature  at  pressure  level  p,  and  Tj(p)  is  given  by 
Eq.  (8c).  Tj(p)  is  the  brightness  temperature  that  the  satellite  would  measure  if  the 
top  of  a  blackbody  cloud  lies  at  pressure  level  p  [which  again  can  be  converted  to  a 
cloud  top  altitude  z  via  the  p(z)  profile]. 

Now  assume  that  the  RTNEPH  has  properly  made  a  cloud  decision  for  a  pixel 
based  on  Eq.  (7).  Consider  the  cloud  top  to  be  at  pressure  level  p.  The  measured 
brightness  temperature  Tsat  for  that  cloudy  pixel  will  be  less  than  the  true  cloud 
temperature  due  to  attenuation  by  the  atmosphere  above  the  cloud.  The  cloud  top 
temperature  Tcld  is  given  by 


Tcld  ~  Tsat  +  ^Tj(pcJd),  (10) 

where  pc]d  is  the  pressure  level  at  the  cloud  top,  and  ATj(pcld)  is  computed  from  Eq. 
(9)  with  p  =  pcid.  Note  that  Eq.  (10)  reduces  to  Eq.  (7)  for  pcid=p9fc. 

The  true  physical  temperature  Tcld  for  a  cloud  top  at  pressure  level  pcjd  is  given 
precisely  by  Eq.  (10)  for  optically  thick  clouds  [that  is,  for  clouds  whose  emissivity  is 
unity;  see  Eq.  (4a)  ].  Equation  (10)  can  be  used  operationally  by  the  RTNEPH  in 
adjusting  cloud  top  temperatures  to  a  more  accurate  value  as  long  as  the  cloud  top 
pressure  is  known.  Obviously,  the  RTNEPH  has  no  way  of  knowing  a  priori  what 
pcId  is.  However,  a  good  estimate  for  pcJd  can  be  found  by  matching  an  uncorrected, 
cloudy  brightness  temperature  Tgat  (from  the  SGDB)  directly  to  a  first-guess  pressure 
Pinit  using  the  upper  air  temperature-height  profile  T(p).  Substituting  p^  for  pcld  in 
Eq.  (10)  yields 


Tcld  ~  Tsat  +  ATj(Pinit)- 
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In  turn,  this  estimate  for  Tcjd  can  be  matched  with  an  improved  estimate  of  pcid, 
again  using  the  T(p)  profile. 


4.3  Improved  Surface  Skin  Temperature  Analysis 

The  computation  of  water  vapor  attenuation  effects  using  radiative  transfer 
theory  has  potential  not  only  for  cloud/no-cloud  decisions  and  improved  cloud  top 
estimates,  but  for  improving  surface  skin  temperature  measurements  as  well. 


For  notational  simplicity,  define  the  normalized  response  function  4>j*(\)  to  be 


6j*(M  = 


<M*) 

Xj2 

f  4>j(k)dX. 

Xji 


where  cj>j(\)  is  the  response  function  for  sensor  j.  Now  consider  the  case  when  the 
RTNEPH  confidently  decides  that  a  pixel  of  brightness  temperature  Tsat  is  cloud  free. 
Using  Eq.  (11),  the  satellite  radiance  as  given  by  Eq.  (4b)  is  written 


Rj  =  /  Rx  VW  dK- 

Xu 


Substituting  (4a)  for  Rx  into  (12)  yields 

Rj  =  f  *,*WB,(T.„)vcd).  +  /  *,*(\)  /  Bk[T(p')](£W  d\.  (13) 

Xj,  P,rc  dp 

Expand  Bx(Tsfc)  [in  the  first  term  on  the  right  side  of  Eq.  (13)  ]  into  a  Taylor 
series  about  the  measured  satellite  brightness  temperature  Tsat  to  obtain 

dBx(Tgat) 

Bk(T.rc)  =  Bk(T„t)  +  (T,[c-T„.)— — -  +  •  ■  ■  . 

For  lTSrc-TSatl  ^~15K,  terms  of  second  order  and  higher  can  be  neglected  since 
their  sum  is  small  in  comparison  to  the  first-order  term.  Substituting  the  truncated 
Taylor  expansion  back  into  Eq.  (13)  yields 

Kr  Kr  3B,(T  t) 

Rj  =  /V(MB,(T.at)TMfcd\  4-  (Tsfc— T8at)  b— d\ 

\j ,  \j,  1,1 


+  JVW/  Bx[T(p')]— 7-dp' d\. 

p.f,  aP 
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Solving  for  Tsfc  finally  yields 


Ri  -  /  <t>j*M 


Tsfc  — 


BJT„tK,1(t+  ;  Bk[T(p')]|r'lP' 

P.fc 


[dk 


Xj2 


+  T, 


JVM  T,„t.-B-T--dx 


sat  > 


dT 


(14) 


where  Tsfc  is  the  desired  surface  skin  temperature,  and  Tsat  is  the  clear-column 
brightness  temperature  measurement.  The  dB/dT  term  in  the  denominator  of  Eq. 
(14)  is  derived  analytically  by  differentiating  Eq.  (2)  with  respect  to  temperature  to 
obtain 


1  qc2k  6  eC2^T 
dT  _  T2  (eC2/xT-l)2 

and  substituting  Tsat  for  T.  The  constants  ^  and  c2  are  as  given  for  Eq.  (2). 

Equation  (14)  can  be  used  in  the  RTNEPH  to  compute  a  surface  skin 
temperature  Tsfc  whenever  a  clear-column  brightness  temperature  Tsat  is  available 
along  with  coincident  upper-air  data  sufficient  to  compute  a  transmittance  profile 

t\(p)- 

Another  method  for  computing  surface  skin  temperatures  from  clear-column 
brightness  temperatures  is  also  applicable  to  RTNEPH  use.  This  method  is  perhaps 
more  straightforward  than  Eq.  (14)  in  obtaining  a  value  for  Tsfc,  but  may  be 
somewhat  more  difficult  to  apply.  Define  the  central  wavelength  for  sensor  j  to  be 
some  \j  that  satisfies 


RXj  =  Rj  for  some  kjX  <  \j  <  \j2,  (16) 

where  Rx.  is  given  by  Eq.  (4a)  with  \  =  kj,  Rj  is  given  by  Eq.  (4b),  and  kjj,  kj2  are  as 
defined  for  Eq.  (3).  The  existence  of  a  kj  is  assured  by  the  Mean  Value  Theorem. 
The  central  wavelength  varies  with  temperature.  The  temperature  dependence  is 
weak  enough,  however,  to  keep  kj  nearly  constant  over  a  typical  LWIR  sensor 
bandwidth  and  typical  terrestrial  temperatures. 

Solving  Eq.  (4a)  with  k  =  kj  for  the  surface  emission  yields 
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> 


(17) 


where  I9at  is  the  measured  clear-column  radiance,  replacing  Rx.  in  Eq.  (4a).  Equation 
(17)  can  be  used  to  calculate  the  surface  emission  whenever  a  cloud-free  radiance 
measurement  is  available  from  satellite.  If  the  satellite  brightness  temperature 
measurement  is  provided  (as  is  the  case  for  the  RTNEPH  SGDB  data),  then  an  Isat 
can  be  computed  using  Eq.  (2)  with  T  —  Tsat  and  \  =  \j. 

Having  computed  Bx.(Tsfc),  it  can  then  be  converted  to  surface  skin  temperature 
Tsfc  in  the  following  way.  Solving  Eq.  (2)  for  T  yields 


which  relates  temperature  T  directly  to  radiance  Bx.  If  Bx  is  replaced  with  Bx.(Tsfc), 
then  the  surface  skin  temperature  Tsfc  is  simply 

c2  _ I 

ci 

X/B, 

where  X.j  is  the  central  wavelength,  and  Bx.  is  computed  using  Eq.  (17)  for  a  cloud-free 
atmosphere. 

The  success  of  Eq.  (18)  in  computing  Tgfc  is  sensitive  to  a  proper  choice  of  the 
central  wavenumber  \j,  which  varies  as  a  function  of  both  the  surface  skin 
temperature  Tsfc  and  the  temperature  profile  T(p).  However,  since  atmospheric 
absorption  is  low,  the  dependence  of  \j  on  Tgfc  is  likely  to  be  stronger  than  that  on 
T(p).  Fortunately  for  LWER  regions,  the  Planck  radiance  varies  slowly  with 
temperature,  so  that  it  is  likely  the  central  wavelength  varies  slowly  as  well.  Perhaps 
a  table  of  central  wavelengths  can  be  constructed  for  just  a  few  surface  skin 
temperature  ranges  with  no  adverse  effect  on  the  accuracy  of  the  results  obtained 
using  Eq.  (18)  (for  example,  \j  =  ll.87p.rn  for  190-225  K;  ll.84p.rn  for  225-275  K; 
and  ll.79p.rn  for  275-330  K).  As  this  report  presents  only  potential  methods  for 
retrieving  improved  skin  temperatures,  the  construction  of  a  table  of  central 
wavelengths  for  the  OLS-T  and  AVHRR  thermal  channels  is  deferred  as  one  of 
several  tasks  of  the  RTNEPH  attenuation  study. 


- p  (18) 

+  1 
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5.  CANDIDATE  INFRARED  TRANSMITTANCE  MODELS  FOR 
WATER  VAPOR  LINE  ABSORPTION 

Having  developed  the  theory  necessary  for  determining  attenuation  corrections  in 
the  RTNEPH,  the  problem  is  reduced  to  computing  atmospheric  transmittances  in 
real-time  using  a  sensible,  accurate,  and  efficient  technique.  Before  proposing 
candidate  models,  it  is  helpful  first  to  discuss  what  the  difficulties  are  in  computing 
reliable  transmittances,  especially  under  real-time  operational  constraints. 

5.1  Narrow  Band  Transmittance  Functions 

Equation  (la)  is  for  monochromatic  radiation  Rx  only.  It  must  be  convolved  with 
the  sensor  response  function  4>j(k),  ^  g'lveQ  by  Eq.  (4b),  before  it  can  be  used  to 
compute  the  radiance  Rj  measured  in  the  spectral  interval  of  a  broadband  LWIR 
instrument  j.  There  is  no  analytic  solution  to  Eqs.  (4a)  and  (4b);  they  must  be 
integrated  numerically. 

To  integrate  Eq.  (4b)  numerically  requires  that  Eq.  (4a)  be  evaluated  for  a  large 
number  of  “closely  spaced”  values  of  wavelength  \  that  fall  between  and  A.j2.  If 
the  separation  A\  between  wavelength  values  has  to  be  very  small  before  accurate 
radiances  can  be  calculated,  this  procedure  could  be  too  cumbersome  for  real-time 
applications.  If  A k  is  too  large,  then  computed  radiance  errors  result  that  are 
comparable  to  or  larger  than  the  noise  in  the  satellite  radiance  measurements 
(Weinreb  and  Hill,  1980). 

Consider  dividing  the  thermal  IR  spectral  band  from  to  Aj2  into  subintervals, 
each  of  which  is  small  enough  so  that  the  clear-column  radiance  Rn  within  the 
interval  n  is  well  approximated  by 


R.  =  Toi,fc  +  /  Bn[T(p')]^dp',  (19a) 

P»lc 

where  the  Planck  radiance  Bn(T)  =  B(\njT),  \n  is  the  central  wavelength  for 
subinterval  n,  and  where 

Tn  =  TIine  ^sb  Tf‘  (19b) 

The  narrow  band  transmittance  Tn  as  given  by  Eq.  (19b)  is  the  average  transmittance 
for  the  spectral  subinterval  n  due  to  water  vapor  absorption.  From  Eq.  (19b)  it  is 
easy  to  see  that  the  narrow  band  water  vapor  transmittance  can  be  treated  in  a 
simple  analytic  fashion  as  a  product  of  the  average  spectral  line  and  the  continuum 
transmittances.  T|ine  is  the  average  transmittance  due  to  water  vapor  line  absorption; 
Tgjj  is  due  to  collisions  among  water  vapor  molecules  (called  the  self-broadened 
continuum);  and  r{  is  due  to  collisions  between  water  vapor  molecules  and  other 
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constituent  molecules  along  the  path  (called  the  foreign-broadened  continuum).  The 
average  transmittance  can  be  expressed  as  the  product  of  Tline,  Tst>,  and  Tf  for  intervals 
A \n  =  \n2  —  Ani  that  are  small  since  the  continuum  transmittances  are  the 
components  of  the  total  transmittance  that  have  a  weak  dependence  on  wavelength. 

By  dividing  up  a  spectral  band  into  several  smaller  subintervals,  Eqs.  (19a)  and 
(19b)  can  be  used  to  determine  the  radiance  Rn  within  each  subinterval  n.  Then  to 
estimate  the  radiance  Rj  measured  over  the  entire  satellite  sensor  j  bandwidth, 
compute  a  weighted  sum 


N 

Rj  =  ^ - ,  (19c) 

2  fn 


where  N  is  the  number  of  subintervals  in  the  domain  from  A.jX  to  \j2  ;  and  where  the 
weights  fn  are  given  by 


/  <MX) dx 

I 


which  is  the  average  of  the  response  function  <f>j  within  the  subinterval  from  knl  to 
\n2.  Equations  (19a)-(19c)  are  approximations  to  Eqs.  (4a)  and  (4b)  that  lend 
themselves  more  easily  to  computation.  The  accuracy  of  Eq.  (19)  is  a  function  of  the 
behavior  of  the  Planck  radiances  B(\n,T)  and  the  transmittance  T(\n)  as 
temperatures  and  wavelengths  change.  Weinreb  and  Hill  (1980)  have  estimated  that 
the  errors  incurred  in  using  the  approximations  (19)  are  no  worse  than  AVHRR 
sensor  noise  ranges.  However,  the  effects  of  the  assumptions  that  went  into 
formulating  Eq.  (19)  will  be  determined  more  carefully  to  ensure  that  no  serious 
errors  are  incurred  when  computing  satellite  radiances  Rj. 

The  monochromatic  atmospheric  transmittance  i\(p)  at  wavelength  \  for  the 
path  between  the  top  of  the  atmosphere  and  any  pressure  level  p  is  given  by 


Tx(p) 


e-Mp)  f 


(20a) 


where  5x(p)  is  the  optical  depth  for  the  atmospheric  path.  For  the  LWIR  region 
where  water  vapor  is  the  primary  absorber,  optical  depth  is  defined  as 


Mp)  =  -- /  kx(p')  q(p')  dP'  » 

g  P 


(20b) 
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where  g  is  the  acceleration  of  gravity,  kx  is  the  cross-section  absorption  coefficient  for 
water  vapor  (area  per  unit  mass),  and  q  is  the  specific  humidity  of  water  vapor 
(dimensionless).  In  general  the  absorption  coefficient  kx  changes  along  an 
atmospheric  path  as  a  function  kx(p,T)  of  pressure  and  temperature  along  that  path. 
The  wavelength  dependence  of  the  absorption  coefficient  is  significant  as  well. 

The  average  line  transmittance  (see  Eq.  (19b)  J  in  the  narrow  band 

\n,  <  \  <  \n2  for  the  path  between  the  top  of  the  atmosphere  and  any  pressure  level 
p  is 


An2 

Tline  =  T7 -  f  TXdk>  W 

Xn  i 

where  rx  is  the  monochromatic  transmittance  due  to  water  vapor  line  absorption. 
Substituting  Eqs.  (20a)  and  (20b)  into  Eq.  (20c)  yields  the  narrow  band  transmittance 
for  the  inhomogeneous  path  between  the  top  of  the  atmosphere  and  any  pressure  level 
P: 


Tline  (p) 


Xn2  0 

—  f  exp  [  kx(p')  q(p')  dp'  ]  d\  . 

aK  Ki  g  P 


(21) 


where  \nl  and  \n2  are  the  endpoints  of  subinterval  n,  and  A\n  =  Xn2— \nl  .  Equation 
(21)  represents  the  only  mathematically  rigorous  way  of  treating  the  dependence  of 
transmittance  on  pressure,  temperature,  and  wavelength.  However  it  is  difficult  to 
evaluate  the  pressure-dependent  integral  in  Eq.  (21)  separately  for  each  wavelength 
inside  the  spectral  subinterval  [Xnl,\n2l  because  of  the  complex  dependence  of 
absorption  coefficient  on  temperature,  pressure,  and  wavelength. 


5.1.1  TRANSMITTANCE  FUNCTION  PARAMETERIZATIONS 

There  are  band  models  available  that  approximate  the  atmospheric  profiles  of 
Tline(p)  ^or  an  inhomogeneous  path  using  transmittance  models  for  homogeneous 
paths  along  which  temperature  and  pressure  are  constant.  Such  models  are  easier  to 
develop  because  the  complicated  (p,T)  dependence  of  the  average  transmittance 
integral  in  Eq.  (21)  is  simplified. 

Consider  a  homogeneous  atmospheric  path  along  which  pressure  p  and 
temperature  T  are  constants  p0  and  T0  ,  respectively.  The  monochromatic  optical 
depth  5X  for  a  wavelength  in  the  LWIR  water  vapor  absorption  region  is  simply  [from 
Eq.  (20b)  ] 


8x(u)  =  kxu> 
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where  the  absorption  coefficient  kx  is  a  function  kx  (po.T0  )  of  wavelength  and  the 
constant  path  pressure  and  temperature,  and  where 


»2 

u  =  J  p(s')ds' 

S| 


is  the  water  vapor  absorber  amount  along  the  layer  path  ds'  from  st  to  s2  .  Then  the 
average  band  transmittance  t  for  the  homogeneous  layer  [from  Eq.  (2 1)]  is  precisely 


Tline  (Po>Tq>u) 


1 

AX 


J  e-kx(p„,T0)u  dx 
kn  1 


(22) 


The  above  integral  has  no  closed  form  and  must  be  evaluated  numerically.  The 
absorption  coefficient  is  an  extremely  complicated,  strongly  varying  function  of 
wavelength.  To  evaluate  Eq.  (22)  accurately  requires  that  the  absorption  coefficient 
kx  be  computed  at  many  wavelengths  between  \nl  and  Xn2  so  that  it  is 
computationally  prohibitive  to  perform  in  a  real-time  or  routine  sense. 

Equations  (19a)  and  (19c)  are  simple  to  compute;  the  biggest  problem  lies  in 
computing  the  average  transmittances  Tn  as  given  by  Eq.  (19b),  and  in  particular  the 
average  line  transmittance  Tline  as  given  by  Eq.  (22).  As  it  stands,  Eq.  (22)  is  no  easier 
to  evaluate  than  Eq.  (4b);  a  computationally  prohibitive  integration  over  some 
bandwidth  is  required  by  both  equations  if  accurate  results  are  to  be  obtained.  The 
only  difference  between  the  two  is  that  in  Eq.  (22),  monochromatic  transmittances 

_  k  u 

e  x  are  averaged,  while  in  Eq.  (4b)  radiances  Rx  are  averaged. 

There  are  models  available  for  computing  reliable  transmittances,  but  they  are 
too  slow  to  be  considered  for  use  by  real-time  applications  such  as  the  RTNEPH.  One 
such  example  is  the  FASCODE  model  developed  at  the  Geophysics  Laboratory. 
FASCODE  is  a  radiative  transfer  program  that  computes  monochromatic 
transmittances  as  given  by  Eq.  (21)  (for  more  FASCODE  information,  see  Clough  et 
al,  1981).  Its  mathematically  accurate  treatment  of  absorption  is  obtained  at  the 
expense  of  routines  that  are  so  computationally  intensive  they  cannot  be  used  in  real¬ 
time.  However,  FASCODE  can  be  used  to  develop,  test,  and  verify  transmittance 
models  for  the  RTNEPH  that  perform  quickly,  but  still  compute  transmittances  to  a 
level  of  accuracy  sufficient  for  the  proposed  cloud  analysis  radiance  techniques  to 
work  successfully.  Two  such  models  are  discussed  in  the  following  paragraphs. 

Multivariate  Polynomial  Regression.  The  water  vapor  line  transmittances  Tline 
can  be  computed  using  a  method  that  requires  less  time  and  memory  than  does  a 
line-by-line  technique  such  as  FASCODE,  but  is  nearly  as  accurate  (Weinreb  and  Hill, 
1980).  The  method  follows  that  of  Weinreb  and  Neuendorffer  (1973),  and  is  a  least- 
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squares  regression  fit  of  the  following  14-parameter  polynomial  to  transmittances  that 
were  calculated  using  a  full-scale,  line-by-'ine  model.  For  a  homogeneous  layer  in 
which  temperature  and  pressure  are  constant,  the  transmittance  can  be  expressed  by 

In  (  -In  7line  )  =  £  C;(\n)  Xi  ,  (23a) 

i=l 

where  Tune  is  the  water  vapor  line  transmittance  averaged  over  subinterval  n,  and 
where 
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(23b) 

where  p,  T,  and  U  are  the  constant  pressure,  constant  temperature,  and  water  vapor 
absorber  amount,  respectively,  along  a  path  in  the  homogeneous  layer.  The 
polynomial  coefficients  Ci(Xn)  are  derived  from  the  least-squares  regression  over  a 
large  sample  of  pre-computed  homogeneous  path  transmittances  using  a  line-by-line 
method  as  specified  by  (22).  Weinreb  and  Hill  (1980)  list  a  set  of  coefficients  in  their 
paper  for  the  3.5  -  4.1  and  10.2  -  12.9pm  regions.  These  spectral  intervals  cover  the 
DMSS  OLS  and  NOAA  AVHRR  thermal  channels. 

Correlated  k  Distribution.  In  the  following  paragraphs  the  correlated  k  method 
is  presented  and  its  strengths  and  limitations  discussed.  The  correlated  k  method 
retains  many  of  the  characteristics  of  an  intensive  line-by-line  transmittance 
calculation,  but  it  computes  transmittances  much  more  quickly  and  efficiently. 
Therefore,  it  is  another  candidate  technique  for  real-time  use  by  the  RTNEPH. 

The  variable  of  integration  in  Eq.  (22)  is  over  the  independent  variable  X,  with 
kx(p0,T0)  being  the  primary  dependent  variable.  Rather  than  integrating  over 
wavelength,  Eq.  (22)  can  be  expressed  as  an  integration  over  values  of  absorption 
coefficients  whereby  the  transmissivity  e  k  is  convoluted  with  the  frequency  of 
occurrence  P(k)  that  absorption  value  k  occurs  within  the  spectral  interval  n,  where 
Xnl  s  X  <  Xn2.  Rewritten  in  this  form,  the  transmittance  [Eq.  (22)  ]  becomes 

00 

Tune  (Po,T0,u)  =  /P(k)e-kudk,  (24) 

o 
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where  P(k)dk  is  the  frequency  of  occurrence  of  absorption  coefficient  between  k  and 
k-fdk  within  the  spectral  subinterval  n.  Equation  (24)  expresses  the  band 
transmittance  as  a  mean  or  “expected”  value  of  transmissivity  between  \n,  and  \n2, 
where  P(k)  represents  a  normalized  absorption  coefficient  density  function  for  the 
band  at  some  pressure  p0  and  temperature  T0.  If  P(k)  is  known  exactly,  then  (24) 
yields  exactly  the  same  result  as  Eq.  (22). 

The  cumulative  absorption  coefficient  distribution  function  g(k)  is  defined  as 

k 

g(k)  =  /P(k')dk',  (25a) 

o 

so  that  g(k)  is  the  fractional  occurrence  of  absorption  coefficient  less  than  or  equal  to 
k.  The  value  of  g(k)  must  lie  between  0  and  1.  Differentiating  Eq.  (25a)  yields 

dg(k)  =  P(k')  dk'  .  (25b) 

Noting  that  when  k=0,  g=0  and  when  k  =  *>,  g  =  1,  Eq.  (24)  can  be  rewritten 
with  g  as  the  variable  of  integration: 


Tline  (PO’^oiU) 


j  e-k.(P..T.)»  dg 


(26) 


where  kg(p0,T0)  is  simply  the  inverse  of  the  function  g(k).  If  kg(p0,T0)  is  known 
exactly,  then  Eq.  (26)  yields  exactly  the  same  result  as  Eq.  (22).  Note  that  Eq.  (26) 
differs  from  Eq.  (22)  only  by  the  variable  of  integration;  g  has  replaced  \. 

The  following  example  illustrates  the  steps  that  transform  Eq.  (22)  into  Eq.  (26). 
Consider  a  simple  spectral  subinterval  n  with  10  equally-spaced  wavelengths  \t,  \2, 
...  ,  \10  that  has  the  simple  absorption  coefficient  spectrum  plotted  in  Figure  2  (top). 
From  Eq.  (22),  the  average  band  transmittance  Tline  for  this  discrete  case  is 


Tline  (Po»T0,u)  =  ~  2  e 


1  -k„Jp0,T0)u 


10m=i 


A\r 


Substituting  the  values  plotted  in  Figure  2  (top)  with  A\m  =  m+1  —  m  =  1  yields 


Tline 


1  r  -k|U  .  -k.,u  -k4u 

—  e  +  e  +  e 
10  1 


e_ksU  +  e_k|U  +  e~k|U 


+  e~k3U  +  e~k8U  +  e-k3U  +  e-k'u  ] 


since  kX|(p0,T0)  =  k!  ,  kX2(p0,T0)  =  k3  ,  ...  ,  kXio(p0,T0)  =  kv  Combining  like 
terms,  the  above  equation  can  be  rewritten  as 
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ABSORPTION  COEFFICIENT  ABSORPTION  COEFFICIENT 


Figure  2.  Pl< 

(Top),  and 


T,i„,  (P..T..U)  -  ^  [  ^"k'“  +  kl"  +  ] 


=  0.4  e  k,u  +  0.4  e  k,u  +  0.1  r  k,u  +  0.1  e  kflU  . 

(27) 

Next  the  average  transmittance  t  will  be  evaluated  using  Eq.  (26),  and  the  result 
compared  to  Eq.  (27).  Before  evaluating  Eq.  (26)  for  the  Figure  2  case,  the  absorption 
spectrum  first  must  be  transformed  from  wavelength  space  to  g  space.  The 
transformation  that  takes  place  between  Eqs.  (22)  and  (26)  is  illustrated  by 
comparing  the  top  and  bottom  plots  of  Figure  2.  What  occurs  as  a  result  of  this 
transformation  is  that  the  absorption  coefficients  are  shuffled  from  their  original, 
random  order  in  wavelength  space  into  an  orderly,  monotonic  arrangement  in  g 
space.  The  smallest  values  of  g  correspond  to  the  smallest  values  of  the  absorption 
coefficient  k,  and  the  largest  values  of  g  correspond  to  the  largest  values  of  k. 

In  order  to  evaluate  Eq.  (26)  for  this  example,  note  from  Figure  2  (bottom)  that 
g,  =  0.1,  g2  =  0.2,  ...  ,  g10  =  1  and  in  general  gm  =  0.1m  .  Therefore  the  average 
band  transmittance  r  using  Eq.  (26)  for  the  simple,  discrete  case  under  consideration 
is 


Tline  (Po,T0,u) 


10 

V 

m=l 


-kKm(Po»Tn)u 

e 


Substituting  in  the  absorption  coefficient  values  for  the  spectrum  plotted  in  Figure  2 
(bottom)  with  Agm  =  gro+1  —  gm  =  0.1(m+l)  -  0.1m  =  0.1  yields 


-  kiu 


-kill 


-k,u 


-k,u 


-k3u  —  kt1u 


t  =  0.1  [  e  +  e  +  e  +  e  +  e  +  e 
+  e  klU  +  e~k3U  +  e-k<u  +  e"^" 


since  kg|(p0,T0)  =  kt  ,  kg2(p0,T0)  =  k,  ,  ...  ,kgin(p0,T0)  =  k8.  Associating  like  terms, 
the  above  equation  can  be  rewritten  as 

*,iDe(Po,T0,u)  =  0.1  [4e~k|U  +  4e_k'’u  4-  e~k,,u  +  e~kflU  ] 

=  0.4e_k'“  +  0.4  e"k’"  +  0.1  e_k’u  +  0.1  e'k-u  , 

(28) 

which  is  precisely  the  transmittance  value  given  by  Eq.  (27).  This  simple  example 
illustrates  that  the  two  forms  [Eqs.  (22)  and  (26)  ]  for  the  transmittance  t  are 
identical. 

The  advantage  of  transforming  the  absorption  spectrum  from  wavelength  space 
to  g  space  lies  in  the  simplicity  of  evaluating  the  average  band  transmittance  integral. 
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In  order  to  accurately  evaluate  Eq.  (22)  (in  wavelength  space),  the  absorption 
coefficient  kx  must  be  computed  at  thousands  of  values  of  A.  within  spectral 
subintervals  \nl  ^  A  s  An2  that  are  only  0.2- 0.4  p.m  in  width.  This  is  because  the 
absorption  coefficient  is  a  quickly  varying  function  of  wavelength,  so  that  sampling  it 
at  relatively  widely-spaced  values  of  A  will  miss  the  fine,  detailed  wavelength 
structure  that  is  so  critical  to  accurate  calculations  of  average  transmittance.  On  the 
other  hand,  Eq.  (26)  (in  g  space)  yields  an  accurate  transmittance  based  upon 
approximately  ten  values  of  g,  because  the  inverse  cumulative  distribution  function 
k(g)  exhibits  a  smooth  monotonic  behavior  that  makes  it  far  less  sensitive  to  sampling 
errors  than  is  k(A). 

This  point  is  illustrated  reasonably  well  by  comparing  the  plots  of  Figure  2. 
There  is  a  more  complex  shape  to  the  top  plot  of  absorption  coefficient  versus 
wavelength,  while  the  bottom  plot  of  absorption  coefficient  versus  g  is  smooth  and 
well  behaved.  The  point  is  even  better  illustrated  by  comparing  the  Figure  3  plots  of 
a  more  physically  plausible  absorption  spectrum  and  its  associated  inverse  cumulative 
absorption  distribution.  Note  how  far  less  critical  it  is  to  sample  the  absorption 
spectrum  in  g  space  (bottom  plot)  as  frequently  as  it  is  in  wavelength  space  (top  plot). 
The  absorption  spectrum  structure  will  likely  be  missed  if  it  is  sampled  at  wavelength 
values  that  are  too  far  apart  from  each  other,  because  the  function  k(\)  changes  so 
quickly  with  small  changes  in  wavelength.  However,  the  inverse  cumulative 
distribution  function  k(g)  is  smooth  and  changes  monotonically  with  g,  so  that  its 
structure  will  be  captured  far  more  effectively  when  it  is  sampled  at  relatively  widely 
spaced  values  of  g.  Thus,  because  fewer  points  are  needed  to  sample  the  absorption 
spectrum,  evaluating  the  average  transmittance  integral  Eq.  (26)  is  computationally 
much  faster  than  having  to  evaluate  Eq.  (22),  while  still  retaining  a  high  accuracy. 

5.1.2  APPLICATION  TO  INHOMOGENEOUS  PATHS 

To  this  point  the  transmittance  models  discussed  have  only  applied  to 
homogeneous  paths.  In  this  section,  the  theory  will  be  presented  for  applying  these 
models  to  inhomogeneous  paths,  that  is,  to  paths  along  which  temperature,  pressure, 
and  water  vapor  absorber  amounts  change. 

Multivariate  Polynomial  Regression.  The  Weinreb  technique  is  readily 
adaptable  to  computing  transmittances  in  the  thermal  OLS  band.  The  technique 
lends  itself  to  real-time,  operational  RTNEPH  applications  and  i3  capable  of  meeting 
all  RTNEPH  water  vapor  attenuation  requirements  [as  specified  by  Eqs.  (19a)-(19c)] 
for  accuracy  and  efficiency.  The  following  subsection  contains  additional  details  on 
how  the  Weinreb  approximation  is  applied  to  inhomogeneous  atmospheres.  However, 
interested  readers  are  strongly  advised  to  refer  to  details  of  the  Weinreb  technique  in 
papers  by  Weinreb  and  Neuendorffer  (1973)  and  Weinreb  and  Hill  (1980). 


-21- 


From  Eq.  (23b)  it  is  easy  to  see  that  the  transmittance  Tune  is  a  function 
TliDe(p,T,U)  of  the  pressure,  temperature,  and  water  vapor  amount  in  some 
homogeneous  layer.  This  is  of  little  practical  value  for  atmospheric  paths  that  are 
inhomogeneous,  that  is,  along  which  pressure,  temperature,  and  mixing  ratios  change 
significantly  between  the  top  of  the  atmosphere  and  the  ground.  Consider  dividing  up 
such  a  path  into  a  sequence  of  homogeneous  layers  each  of  which  has  p,  T,  and  U 
constant.  If  the  transmittances  Tiine  were  not  a  band  average  but  rather  truly 
monochromatic  for  some  wavelength  X,  then  the  total  atmospheric  transmittance 
over  the  inhomogeneous  path  from  the  surface  to  the  top  of  the  atmosphere  would  be 
the  product 


Tline(^) 


II  Tline  (  P/>  T„  U,,  X  ) 
1=1 


where  Tline(X)  is  the  monochromatic  transmittance  due  to  water  vapor  line 
absorption,  and  where  P;  ,T/ ,  and  U;  are  the  homogeneous  pressure,  temperature,  and 
water  vapor  amount,  respectively,  for  layer  /;  and  where  L  is  the  number  of  layers 
into  which  the  atmosphere  has  been  subdivided. 

If  Eq.  (23)  were  used  to  compute  the  layer  transmittances  Tline(p/ ,T;  ,U/  ,X)  on 
the  right  side  of  the  above  equation,  the  total  atmospheric  attenuation  Tline(X)  would 
be  miscalculated.  This  is  because  the  transmittance  along  the  inhomogeneous  path  is 
simultaneously  dependent  on  all  the  conditions  (p^T^Uj),  (P2>T2,U2),  ...  , 

(pL,TL,UL).  The  technique  that  is  used  for  inhomogeneous  paths  is  as  follows. 

Equation  (23)  establishes  a  relation 

Time  =  P,  T,  U  )  (29a) 

between  transmittance  and  pressure,  temperature,  and  water  vapor  amount.  There 
is  also  a  relation 


u  =  U(  p,  T  )  (29b) 

that  equates  average  transmittance  Tiine  to  water  vapor  amount  U  for  a  given  p  and 
T.  The  relation  (29b)  is  determined  simply  by  solving  Eq.  (23)  for  U.  With  these 
functional  relationships  established,  the  Weinreb  approximation  to  inhomogeneous 
paths  can  be  described. 

Consider  an  atmosphere  composed  of  L  homogeneous  layers  each  of  some 
pressure,  temperature,  and  absorber  amount  p / ,  T / ,  and  U/,  respectively.  The  total 
transmittance  from  the  top  of  the  atmosphere  through  the  first  layer  (the  top  layer)  is 
simply 
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Tline(  Pl>  Tlt  Ux  ), 


(30a) 


which  for  notational  simplicity  shall  be  called  tx.  For  the  transmittance  t2  through 
both  the  first  and  the  second  layers,  there  is  a  bit  of  a  dilemma.  Equation  (23)  only 
computes  transmittances  for  single  values  for  p,  T,  and  U;  what  is  needed  now  is  the 
transmittance  for  a  path  along  which  p,  T,  and  U  change. 

The  total  amount  of  water  vapor  in  the  path  through  layers  1  and  2  is  given  by 
the  sum  Ux  +  U2-  If  the  pressure  and  temperature  through  layers  1  and  2  are 
approximated  by  p2  and  T2,  then  the  total  transmittance  t2  for  the  path  through  the 
first  two  layers  generally  will  be  underestimated  if  it  is  computed  as 

t2  =  ^line(  P2>  T2.  Ux+U2  )•  (30b) 

However,  there  is  a  value 

U2*  =  U2  +  AU2 

with  which  the  total  transmittance  t2  under  both  conditions  px,  Tx,  and  p2,T2  can  be 
well  approximated  as  a  function  t2  =  Tj|ne(  p2,  T2,  U2*).  U2*  is  called  the  scaled 
water  vapor  amount,  and  is  the  sum  of  the  known  amount  of  water  vapor  U2  in  layer 
2  and  the  unknown  AU2  whose  value  must  be  determined.  For  the  path  between  the 
top  of  the  atmosphere  and  the  bottom  of  layer  2, 


T2  =  Tline(  P2)  T2,  U2+AU2  ). 

But  at  the  boundary  between  layer  1  and  layer  2  the  total  atmospheric  transmittance 
is  Tx.  Also,  at  the  very  top  of  layer  2,  the  water  vapor  amount  U2  is  zero,  so  that  the 
above  equation  is  simply 


Ti  =  Tline(  P2,  T2,  AU2).  (31a) 

The  transmittance  tx  is  given  by  Eq.  (30a);  thus  the  only  unknown  in  the  above 
relation  is  the  value  AU2.  Using  Eq.  (29b)  it  is  simply 

AU2  =  U(t1,P2,T2),  (31b) 

which  states  that  AU2  is  precisely  the  amount  of  water  vapor  that  specifies  the 
transmittance  tx  through  the  first  layer,  except  under  the  conditions  p2,  T2  of  the 
second  layer.  Again,  the  scaled  absorber  amount  U2*  is  the  sum  U2  +  AU2  of  the 
water  vapor  amount  in  layer  2  and  the  scaled  amount  for  layer  2.  Finally,  the  total 
transmittance  t2  through  layers  1  and  2  is  given  by 
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t2  —  TIme(  P2»  T2,  U2+AU2  ). 


(32) 


In  general  Uj  ^  AU2,  indicating  that  the  total  transmittance  t2  between  layers  1  and 
2  is  underestimated  when  using  Eq.  (30b)  as  opposed  to  Eq.  (32). 

For  an  atmosphere  that  has  been  subdivided  into  L  homogeneous  layers,  a 
recurrence  relation  can  be  set  up  whereby 


t /  =  Tline(  p,,  T/(  U/+ AU/)  ,  and 

(33) 

AU/  =  U(  T,_lf  p/(  T,)  for  /  =  2,  3,  L, 

where  I-!  =  Tune(  pj  ,  Tj  ,  Uj  )  as  given  by  Eq.  (23).  The  above  procedure  is  repeated 
L  times  to  arrive  at  the  average  transmittance  =  tl  in  subinterval  n,  as  required 
by  Eq.  (19b),  where  t i  is  the  average  transmittance  due  to  water  vapor  line 
absorption  for  the  inhomogeneous  path  between  the  top  of  the  atmosphere  and  the 
ground. 

Correlated  k  Distribution.  Equation  (26)  can  be  accurately  evaluated  for  a 
homogeneous  path  along  which  pressure  and  temperature  are  constant.  However, 
pressure  and  temperature  variations  strongly  influence  the  shape  of  absorption 
spectra  and  hence  the  associated  cumulative  absorption  distribution.  Thus  k(g)  is  a 
function  k(g,p,T)  of  pressure  and  temperature.  In  order  to  calculate  transmittance  for 
a  path  along  which  pressure  and  temperature  are  not  constant,  the  following 
correlated  absorption  method  is  used. 

Consider  a  spectral  subinterval  [Xni,\n2]  within  which  the  absorption  lines  are 
equally  spaced  and  of  equal  intensity  as  shown  in  Figure  4.  Consider  also  an 
atmosphere  composed  of  L  homogeneous  layers  within  which  pressure  p/  and 
temperature  T t  are  constant.  The  difference  between  the  absorption  spectra  from  one 
layer  to  the  next  is  due  only  to  differences  in  layer  pressure  and  temperature.  For 
notational  simplicity,  let 


k/(g)  =  k  (  P/,  T,,  g  ) 

for  the  Ith  layer,  where  /  =  1,  2,  ...,  L.  (In  this  notation  the  subscript  /  denotes  the 
pressure-temperature  dependence  of  the  absorption  coefficient.) 

For  this  idealized  case,  there  is  some  wavelength  K  which  is  associated  with  a 
value  k^1)  of  the  absorption  coefficient  in  the  first  (topmost)  layer,  there  is  always  the 
same  set  of  precisely  known  values  k^2\  k^3^,  ...,  k^  of  absorption  coefficients  in 
layers  2,  3,  ...,  L,  respectively,  that  are  associated  with  \.  Thus  there  is  a  one-to-one 
mapping  and  a  complete  correlation  between  the  absorption  coefficients  (hence  the 
name  “correlated  k”)  in  the  first  layer  and  those  in  successively  lower  layers  (see 
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ABSORPTION  COEFFICIENT  k(X,p,T) 


X„,  WAVELENGTH  Xn2 


Figure  4.  Plot  of  an  Idealized  Absorption  Spectrum  Profile  Illustrating 
the  Conditions  Necessary  for  the  Correlated  k  Transmittance  Model  to 
be  Valid  for  an  Inhomogeneous  Atmospheric  Path 
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Figure  4).  It  follows  that  the  fraction  of  wavelengths  with  absorption  coefficient 
values  less  than  in  layer  1  must  be  the  same  as  the  fraction  of  wavelengths  with 
absorption  coefficients  less  than  k^2',  k^,  . . .  ,  k^  in  layers  2,  3,  L,  respectively. 
Mathematically,  this  is  expressed  as 

gi  (  k^  )  =  g2  (  k(2)  )  =  ...  =  gL(  k(i>  )  .  (34) 

Simply  stated,  the  percentile  rank  of  absorption  coefficient  at  a  particular  wavelength 
within  a  spectral  subinterval  remains  unchanged  from  one  layer  to  che  next;  that  is, 
absorption  coefficient  rank  is  independent  of  pressure  and  temperature.  For  example, 
if  absorption  coefficient  value  k'1)  is  greater  than  80  percent  of  the  k  values  in  layer  1, 
then  k^  is  greater  than  80  percent  of  the  k  values  in  layer  /,  for  /  =  2,  3,  ...,  L.  Under 
these  assumptions  an  inverse  cumulative  distribution  function  can  be  constructed  for 
the  optical  depth  of  an  inhomogeneous  path  of  L  layers  by  simply  adding  the  optical 
depths  for  each  homogeneous  layer  as  a  function  of  g;  it  is  given  by 

8(g)  =  2  k/(g)u/>  (35a) 

l—i 


where  k/(g)  is  the  inverse  of  the  cumulative  absorption  coefficient  distribution 
function  g,(k)  for  layer  l,  and  the  absorber  amount  u/  for  layer  /  is  [from  Eq.  (20b)] 

u /  =  —  q(p/)  Ap , ,  (35b) 

9 

where  g  is  the  acceleration  due  to  gravity  (the  italic  g  denotes  the  acceleration 
constant  as  opposed  to  the  correlated  k  transformation  variable  “g”),  and  where  Apj 
is  the  thickness  for  layer  /  (in  pressure  units,  positive-definite). 

For  this  idealized  case,  then,  the  average  transmittance  t  for  the  inhomogeneous 
path  between  the  top  of  the  atmosphere  and  some  pressure  level  p L  is  [from  Eq.  (26)] 

1  L 

Tlme(pJ  =  J  exP  [  -  2  k/(g)  u/  ]  dg  •  (35c) 

0  1=1 

Correlated  k  methods  work  best  within  absorption  bands  that  satisfy  the 
idealized  conditions  illustrated  in  Figure  4,  and  stated  by  Eq.  (34).  The  extent  to 
which  Eq.  (34)  is  valid  dictates  the  accuracy  of  Eq.  (35c).  In  other  words,  how  well 
the  absorption  coefficients  k(g)  are  correlated  from  one  atmospheric  layer  to  the  next 
is  the  critical  assumption  that  determines  in  general  whether  an  accurate  average 
transmittance  value  can  be  computed  using  a  correlated  k  technique. 

In  the  LWIR  water  vapor  absorption  spectrum,  line  strengths  and  spacings  are 
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not  uniform  as  they  are  for  the  idealized  case  in  Figure  4.  Temperature  affects  the 
strengths  of  some  lines  more  strongly  than  others,  especially  for  lines  near  the  wings 
of  molecular  absorption  bands.  Hence,  line  strengths  vary.  In  addition,  the  distance 
between  lines  is  not  even,  causing  lines  to  overlap  in  non-uniform  ways  as  pressure 
and  temperature  change.  For  a  value  of  absorption  coefficient  in  layer  1  there 
exists  not  one  but  rather  a  range  of  values  of  absorption  coefficients  associated  with 
k^  in  any  subsequent  layer.  Because  of  these  effects  the  absorption  coefficients  are 
not  exactly  correlated,  as  Eq.  (34)  specifies  they  should  be. 

However,  the  degiee  of  lack  of  correlation  is  small;  that  is,  the  range  of  k  values 
associated  with  k^  is  small  so  that  Eq.  (34)  remains  a  good  approximation.  The  high 
degree  of  correlation  among  absorption  coefficients  from  one  atmospheric  layer  to  the 
next  has  allowed  the  correlated  k  method  to  be  successfully  applied  to  other 
problems. 

From  the  previous  subsection,  recall  that  for  a  homogeneous  layer  the  advantage 
of  transforming  the  original  absorption  spectrum  from  wavelength  space  to  g  space 
lies  in  the  simplicity  of  the  spectrum  shape  (refer  back  to  Figure  3).  In  g  space  the 
absorption  spectrum  k(g)  is  smoother  and  hence  is  less  sensitive  to  discrete  sampling 
errors  than  it  is  in  wavelength  space.  By  dividing  up  g  space  into  several  smaller 
subdomains,  Eq.  (35c)  can  be  written  in  a  quadrature  form  as 

Tline(Pz)  =  2  exP  (  “  S  k/(&)  u/  1  Ag5 ,  (36) 

i=0  1—1 

where  I  is  the  number  of  times  the  absorption  spectrum  k;(gj)  is  to  be  sampled  within 
the  interval  Osg<l,  Equation  (36)  can  be  used  to  compute  the  average  water 
vapor  transmittance  [see  Eq.  (19b)]  in  the  spectral  subinterval  Xnl  <  X  <  Xn2  for 
the  inhomogeneous  atmospheric  path  between  the  top  of  the  atmosphere  and  any 
pressure  level  pL. 

Techniques  have  been  presented  to  compute  narrow  band  water  vapor  line 
transmittances  Tune  [from  Eq.  (19b)]  that  are  fast  and  accurate.  The  next  section  is 
devoted  to  presenting  techniques  that  will  be  used  to  calculate  the  continuum 
transmittances  Tsb  and  Tf  [also  from  Eq.  (19b)  ]. 


6.  INFRARED  TRANSMITTANCE  MODELS  FOR  WATER  VAPOR 
CONTINUUM  ABSORPTION 

Typically  two  types  of  water  vapor  continuum  absorption  are  treated:  self- 
broadened,  which  accounts  for  attenuation  due  to  collisions  between  water  vapor 
molecules,  and  foreign-broadened,  which  accounts  for  collisions  among  water  vapor 
molecules  and  other  atmospheric  constituents.  The  self-broadened  absorption  is 
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dependent  upon  the  partial  pressure  of  water  vapor,  while  foreign-broadened  is 
dependent  upon  the  partial  pressure  of  the  dry  atmosphere  (Weinreb  and  Hill,  1980). 

6.1  Self-Broadened  Continuum 

The  self-broadened  average  transmittance  Tsb  in  spectral  subinterval  n  for  a  path 
between  the  top  of  the  atmosphere  and  some  pressure  level  p  is  given  by  (Weinreb 
and  Hill,  1980) 


-In  [  Tgb(p,\n)  ]  = 


P  T0(— 

5.41  X  1013  C°(Xn)  sec0  J*  p' r2  e  T 

o 


298  J 


dp', 


(37) 


where  p  is  the  atmospheric  pressure  that  specifies  the  end  of  the  path;  Xn  is  the 
central  wavelength  for  the  spectral  subinterval  n  from  Xnl  to  Xn2;  C°(Xn)  are 
coefficients  determined  for  each  subinterval  n;  0  is  the  angle  between  the  atmospheric 
path  and  the  local  vertical  (the  zenith  angle);  r  is  the  mass  mixing  ratio  of  water 
vapor  (  g/kg,  dimensionless  )  at  pressure  level  p';  T0  is  1800  K  for  the  10.2  —  12.9p.rn 
region  and  1300  K  for  the  3.5  —  4.1p.m  region;  and  T  is  the  atmospheric  temperature 
at  pressure  p'. 

Equation  (37)  can  be  written  in  quadrature  form  for  a  sequence  of  L  homogeneous 
layers  in  spectral  subinterval  n  as 


L  To  (  _  x  ) 

-In  l  T9b(PL  Aj  ]  =  5.41X  1013  C°(Xn)  sec0  P/r2e  T'  298  Ap, , 

/=  l 

(38) 

where  P/,r/,  and  T/  are  the  constant  pressure,  mixing  ratio,  and  constant 
temperature,  respectively,  for  homogeneous  layer  /;  Ap;  is  the  thickness  (in  pressure 
coordinates)  of  layer  /;  and  where  all  other  variables  are  as  in  Eq.  (37).  The 
coefficients  C°(Xn)  are  listed  in  Weinreb  and  Hill  (1980)  for  each  subinterval  n. 

6.2  Foreign-Broadened  Continuum 

The  foreign-broadened  water  vapor  continuum  average  transmittance  Tf  in 
spectral  subinterval  n  for  a  path  between  the  top  of  the  atmosphere  and  some 
pressure  level  p  is  given  by  (Weinreb  and  Hill,  1980) 


p 

— In  [  Tf  (p,\n)  ]  =  4.04  x  1015  C°(Xn)  sec0  f  p'  r  dp',  (39) 

o 

where  p  is  the  atmospheric  pressure  that  specifies  the  end  of  the  path;  Xn  is  the 
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central  wavelength  for  the  spectral  subinterval  n  from  \nl  to  ka2  [the  same  as  in  Eq. 
(37)];  C°(\n)  are  coefficients  determined  for  each  subinterval  n;  0  is  the  angle  between 
the  atmospheric  path  and  the  local  vertical  (the  zenith  angle);  and  r  is  the  mass 
mixing  ratio  of  water  vapor  (  g/kg,  dimensionless  )  at  pressure  p'.  In  longwave 
thermal  IR  window  regions,  the  foreign-broadened  continuum  absorption  is  negligible 
in  comparison  with  self-broadened  and  line  absorptions;  thus  Eq.  (39)  is  ignored  for 
10.2  -  12.9p.m. 

Equation  (39)  can  be  written  in  quadrature  form  for  a  sequence  of  L  homogeneous 
layers  in  spectral  subinterval  n  as 


“In  [  Tf  (Pz,  ,ka)  ]  =  4.04X1015  C°(kn)  secG  2  P/r,Ap;, 

/=i 

(40) 

where  p/,r/,  and  are  the  constant  pressure,  mixing  ratio,  and  constant 
temperature,  respectively,  for  homogeneous  layer  /;  Ap/  is  the  thickness  (in  pressure 
coordinates)  of  layer  /;  and  where  all  other  variables  are  as  in  Eq.  (39).  Remember, 
Eq.  (40)  is  for  the  3.5  -  4.1pm  spectral  interval  only;  foreign-broadened  continuum 
absorption  is  ignored  in  the  10.2  -  12.9pm  region,  but  it  cannot  be  ignored  in  the 
3.5  -  4.1pm  region  (Weinreb  and  Hill,  1980). 

Methods  have  been  discussed  that  allow  for  quick,  accurate  computations  of  the 
average  water  vapor  transmittance  Tn  [as  given  by  Eq.  (19b)]  over  spectral 
subinterval  n  (from  \nl  to  \n2)  for  an  inhomogeneous  atmospheric  path.  Using  Eq. 
(19b),  Tn  is  well  approximated  as  the  product  of  the  water  vapor  line  transmittance 
Tline  as  given  by  either  Eq.  (33)  (for  the  multivariate  regression  technique  of  Weinreb 
and  Hill)  or  Eq.  (36)  (for  the  correlated  k  method),  the  self-broadened  continuum 
transmittance  Tsb  as  given  by  Eq.  (38),  and  the  foreign-broadened  continuum 
transmittance  Tf  as  given  by  Eq.  (40). 


7.  OVERLAPPING  ABSORPTION  BANDS 

The  narrow  band  transmittance  Tn  as  given  by  (19b)  is  affected  by  absorption  due 
to  constituents  other  than  water  vapor,  such  as  molecular  nitrogen  and  the  uniformly 
mixed  gases.  In  formulating  (19b)  it  was  assumed  that  water  vapor  absorption  is 
dominant  in  LWIR  window  regions.  This  is  certainly  the  case  for  the  OLS  and 
AVHRR  thermal  IR  channels;  however,  the  absorption  bands  of  molecular  nitrogen 
and  other  uniformly  mixed  gases  overlap  the  water  vapor  absorption  bands  within 
these  channels  and  have  small  effects  on  total  transmittance. 

The  transmittance  for  a  particular  wavelength  k  is  affected  by  the  quantity 
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and  type  of  absorbers  present  in  an  atmospheric  path.  Assume  for  illustration  that 
there  are  three  types  of  absorbers  a,  b,  and  c  in  a  path  along  which  transmittance  is 
to  be  computed.  The  total  monochromatic  transmittance  along  that  path  is  given  by 

=  Ta(A)  Tb(A)  Tc(\), 

where  xa(A)  is  the  transmittance  due  to  absorber  type  a,  etc.  Inserting  this  form  into 
Eq.  (20c)  yields 


T 


n 


X 


n  2 


/  Ta(\)  Tb(\)  Tc(\)  dA 

^ii  i 


A  A 


n 


Thus  the  average  total  transmittance  in  spectral  band  n  is  the  average  of  the 
products  of  the  monochromatic  transmittances  for  each  absorber  type  along  the  path. 
If  the  spectral  band  is  chosen  to  be  small  enough  then  the  average  total  transmittance 
can  be  well  approximated  by 


Tn  ~ 


_ _  i 


f  Ta(A)  dA  /  Tb(A)  dA  /  tc(A)  dA 


An  i 


AA, 


In  general  the  above  equation  is  not  precisely  true  except  for  monochromatic 
transmittances  (that  is,  when  Anl  =  An2)  or,  for  A^^A^,  if  the  constituents  a,  b,  and 
c  have  absorption  characteristics  that  are  completely  uncorrelated  with  each  other. 
However,  the  product  of  the  average  individual  transmittances  over  the  interval 
generally  approximates  well  the  average  product  of  the  monochromatic 
transmittances. 

Thus,  assuming  that  the  spectral  interval  n  is  properly  chosen,  the  average  total 
transmittance  Tn  as  given  by  Eq.  (19b)  can  be  expressed  as  the  product  of  the  average 
transmittances  for  the  individual  absorbers.  From  the  above  equation,  then,  the 
average  transmittance  for  spectral  band  n  is 

tH20  tN  Tgas  >  (41) 


where 


tH20  =  Tline  Tsb  Tf 

is  the  average  transmittance  for  spectral  band  n  due  to  water  vapor  line  absorption 
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[see  Eq.  (19b)];  ts5  is  due  to  collisions  among  water  vapor  molecules  (the  self- 
broadened  continuum);  Tf  is  due  to  collisions  between  water  vapor  molecules  and 
other  constituent  molecules  along  the  path  (the  foreign-broadened  continuum);  tn  is 
due  to  collision-induced  molecular  nitrogen  absorption;  and  Tgag  is  the  average  band 
transmittance  due  to  uniformly  mixed  gases  (carbon  dioxide,  nitrous  oxide,  carbon 
monoxide,  methane,  and  molecular  oxygen). 

There  are  techniques  similar  to  Eqs.  (38)  and  (40)  that  are  available  for  use  in 
computing  tn  and  Tgas ;  examples  are  presented  in  Weinreb  and  Hill  (1980).  As 
previously  mentioned,  water  vapor  is  the  primary  absorber  for  the  thermal  infrared 
wavelengths.  As  far  as  the  RTNEPH  is  concerned,  attenuation  by  molecular  nitrogen 
is  significant  only  in  the  region  near  4.3p.m,  and  the  uniformly  mixed  gases  absorb 
weakly  in  both  the  3.7  and  10-12jjim  regions.  Thus  for  RTNEPH  purposes  these 
effects  will  likely  be  neglected,  since  they  are  typically  small  in  comparison  to 


8.  DISCUSSION  AND  PLANS 

Techniques  have  been  presented  that  can  be  used  to  compute  the  expected  IR 
water  vapor  attenuation  along  a  given  atmospheric  path.  Such  computations  are 
necessary  for  cloud/no-cloud  threshold  techniques  (such  as  the  RTNEPH  IR 
processor)  that  compare  a  measured  IR  brightness  temperature,  which  includes  the 
effects  of  atmospheric  water  vapor  attenuation,  with  a  corresponding  surface  skin 
temperature.  In  this  section,  study  plans  are  outlined  to  show  how  these  attenuation 
models  will  be  incorporated  into  the  RTNEPH. 

Before  the  homogeneous  path  transmittance  model  of  Weinreb  and  Hill  (see 
Section  5.1)  can  be  incorporated  into  the  RTNEPH,  several  minor  efforts  must  first 
be  completed.  Recall  from  Eqs.  (19a)-(19c)  that  the  basic  method  for  computing  an 
expected  clear-column  IR  radiance  is  initially  to  subdivide  the  thermal  band  of 
interest  into  smaller  subintervals,  compute  an  average  radiance  for  each  subinterval 
[Eqs.  (19a)  and  (19b)],  and  then  perform  a  weighted  average  of  the  subinterval 
radiances  to  get  the  total  band  radiance  [Eq.  (19c)].  The  band  radiance  is  then 
converted  to  an  expected  clear-column  brightness  temperature  using  Eq.  (18). 
Currently  the  homogeneous  path  technique,  originally  developed  for  the  relatively 
narrow  AVHRR  thermal  channels,  divides  the  broadband  DMSS  OLS-Thermal  band 
into  eight  subintervals.  A  sensitivity  study  must  be  performed  to  ensure  that  these 
eight  subintervals  will  allow  for  OLS-T  radiance  calculations  of  sufficient  accuracy. 
(The  complexity  of  the  IR  sensor  response  function  along  with  the  spectral  variation 
of  the  Planck  function  across  the  band  determines  how  many  subintervals  need  to  be 
chosen.)  For  purposes  of  comparison,  errors  in  clear-column  brightness  temperature 
of  less  than  0.5  K  from  the  line-by-line  FASCODE  brightness  temperature  will  be  the 
accuracy  benchmark  since  0.5  K  will  be  roughly  the  thermal  resolution  of  an  8-bit  IR 
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SGDB. 

The  Weinreb  and  Hill  continuum  calculations  [Eqs.  (38)  and  (40)]  will  be  replaced 
by  the  latest  FASCODE  water  vapor  continuum  model.  This  task  essentially  involves 
updating  the  coefficients  C°(Xn)  that  appear  in  Eqs.  (30)  and  (32),  and  formulating 
the  (p,T)  effects  slightly  differently  (see  Clough  et  al,  1981).  Also,  the  water  vapor 
line  absorption  regression  coefficients  Cj(\n)  in  Eq.  (23)  may  be  recomputed  (via 
multivariate  regression  to  FASCODE  water  vapor  line  absorption  transmittances)  if 
it  is  determined  that  the  line  transmittances  do  not  match  those  of  FASCODE  to 
within  1  or  2  percent.  Such  a  tight  constraint  is  necessary  because  small  errors  in 
transmittance  can  manifest  themselves  as  larger  errors  in  computed  brightness 
temperature  via  the  differential  transmission  term  drn  /  dp  [see  Eq.  (19a)]. 

With  these  issues  addressed,  attention  can  then  be  focus  on  the  implementation 
of  a  candidate  transmittance  model  (either  the  multivariate  polynomial  regression  or 
correlated  k)  into  the  RTNEPH  operational  framework.  Validation  runs  will  test  the 
accuracy  of  the  transmittance  model  against  ground  truth  measurements  of  dear- 
column  brightness  temperature  for  a  known  atmospheric  path.  Among  the  data 
available  for  such  tests  include  maritime  radiosonde  data  in  conjunction  with  a 
measured  sea  surface  temperature  and  coincident  IR  brightness  temperature 
measurements.  (Such  tests  cannot  be  performed  over  land  because  of  the 
unavailability  of  routine  ground  truth  skin  temperature  measurements.)  Other 
verification  test  data  sources  include  coincident  measurements  of  surface  skin 
temperature,  upper  air  temperature  and  water  vapor  profiles,  and  satellite-based  IR 
radiance  measurements  obtained  during  the  FIFE  ’89  summer  land- 
surface/atmosphere  field  experiment  over  the  Kansas  plains. 

Checks  will  also  be  performed  against  FASCODE  runs  to  ensure  the 
transmittance  models  are  accurate.  Initial  comparisons  have  checked  out  to  well 
within  the  0.5  K  tolerance,  but  several  more  FASCODE  runs  will  be  made  over  a 
range  of  globally-expected  atmospheric  paths,  ranging  from  very  dry  at  low  viewing 
angles  to  very  moist  at  high  viewing  angles.  In  terms  of  accuracy  it  is  expected  that 
the  attenuation  model  will  be  more  than  adequate  for  the  RTNEPH.  It  is  more  likely 
that  any  observed  shortcomings  in  the  transmittance  calculations  will  be  due  to 
’  inaccurate  and/or  insufficient  surface  skin  temperatures  and  atmospheric  profile  data 

•  available  to  the  model. 

Once  the  attenuation  model  is  ready,  the  first  step  will  be  to  compare  its  water 
vapor  attenuation  correction  calculations  to  those  currently  performed  by  the 
RTNEPH.  A  study  is  now  underway  that  computes  the  water  vapor  attenuation  for 
upper  air  stations  throughout  North  America  over  a  period  of  several  months.  The 
chosen  radiosonde  stations  are  plotted  in  Figure  5.  Note  that  a  wide  variety  of 
tropical,  maritime,  arctic,  and  continental  stations  have  been  selected.  The  plan  is  to 
collect  radiosonde  data  from  these  stations  and  run  three  LWIR  water  vapor 
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Figure  5.  Radiosonde  Stations  Selected  for  the  Water  Vapor  Attenuation  Study 


attenuation  codes  on  their  atmospheric  profile  data:  the  candidate  RTNEPH 
technique  (multivariate  polynomial  regression  or  correlated  k),  LOWTRAN  (a 
parameterized  version  of  FASCODE),  and  the  current  operational  RTNEPH  “lookup 
table”  scheme.  The  current  and  proposed  RTNEPH  attenuation  correction 
techniques  will  then  be  directly  compared.  Through  this  study  the  following  issues 
will  be  directly  addressed: 

1.  From  a  computational  point  of  view,  how  complex  does  the 
transmittance  model  need  to  be  in  order  to  achieve  a  computed 
equivalent  radiance  error  less  than  0.5  K?  Is  the  current  RTNEPH 
attenuation  scheme  accurate  enough  or  is  a  more  sophisticated  scheme 
warranted? 

2.  Is  there  an  observed  and  significant  fine-scale  horizontal  variation  in 
water  vapor  attenuation  that  is  captured  by  the  resolution  of  the  data 
in  this  study  but  that  cannot  be  captured  by  the  RTNEPH  HIRAS 
upper  air  database?  If  this  is  so,  does  the  whole-mesh  upper  air 
temperature  profile  somehow  need  to  be  updated  throughout  the 
boundary  layer  at  eighth-mesh  resolution  using  the  local  surface  skin 
temperature  measurement? 

3.  Is  there  a  significant  loss  of  information  if  significant  level  data  are  left 
out  of  the  RTNEPH  upper  air  temperature  profiles  (that  is,  if  only 
mandatory  level  data  are  used)? 

A  preliminary  estimate  can  likely  be  made  through  this  study  of  how  the  proposed 
transmittance  computations  might  affect  the  RTNEPH  IR  processor  cloud/no-cloud 
analysis. 

Currently,  the  RTNEPH  uses  four  lookup  tables  to  help  determine  water  vapor 
attenuation  corrections  as  a  function  of  1)  IR  grayshade,  2)  sensor  viewing  angle,  3) 
whether  the  point  of  interest  lies  within  the  tropics,  and  4)  day/night,  sea/land,  and 
geographic  location  (that  is,  RTNEPH  box).  These  complicated  corrections  are 
necessary  because  of  biases  in  the  spectral  response  for  a  particular  IR  sensor  as  it 
ages  onboard  a  spacecraft,  data  truncation  and  viewing  angle  correction  procedures 
that  the  IR  SGDB  software  invokes  on  the  raw  satellite  data  before  it  is  processed  by 
the  RTNEPH,  and  because  OLS-Thermal  data  are  not  precisely  calibrated.  The 
attenuation  correction  techniques  proposed  in  this  paper  will  eliminate  the  need  for 
the  first  three  RTNEPH  lookup  table  correction  factors,  replacing  them  with  a 
physically  rigorous  and  robust  technique. 

In  and  of  itself,  it  is  unrealistic  to  expect  that  the  insertion  of  a  physically  sensible 
method  of  computing  water  vapor  attenuation  correction  factors  will  eliminate  all  of 
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the  errors  in  the  RTNEPH  cloud/no-cloud  analysis.  The  attenuation  computations 
may  very  well  be  accurate  for  the  atmospheric  profile  data  that  were  used  to  compute 
them.  However,  if  these  atmospheric  data  are  themselves  inaccurate,  errors  in  the 
cloud  analysis  will  likely  occur.  A  need  for  including  some  type  of  bias  correction  will 
still  exist  even  with  the  new  attenuation  model.  Thus  some  type  of  lookup  table  that 
is  a  function  of  day/night,  sea/land,  and  geographic  location  will  likely  be  retained 
because  of  biases  in  the  input  data  and  uncertainties  in  sensor  degradation  effects. 

The  surface  radiance  term  [first  term  on  the  right  side  of  Eq.  (19a)]  is  dominant 
in  most  IR  upwelling  radiance  calculations.  The  surface  radiance  is  strongly 
dependent  on  the  surface  skin  temperature.  Therefore,  accurate  values  of  surface 
skin  temperature  are  important  in  computing  clear-column  radiances.  To  this  end, 
improvements  are  being  proposed  and  developed  for  the  surface  skin  temperature 
model  used  by  the  RTNEPH  IR  processor  threshold  technique.  Presumably  the 
characteristics  of  the  new  skin  temperatures  will  prompt  a  significant  and  nontrivial 
retuning  effort  of  the  current  RTNEPH  water  vapor  attenuation  correction  factors. 
At  the  same  time  GL  is  planning  to  develop  and  implement  a  new  method  for 
calculating  these  correction  factors.  If  one  upgrade  (for  either  temperature  or 
attenuation)  becomes  available  and  is  implemented  before  the  other,  then  a  two-step 
process  is  envisioned  whereby  two  long-term  retuning  periods  will  result;  one  after 
implementation  of  the  new  surface  skin  temperature  model  and  another  after 
implementation  of  the  new  attenuation  correction  model  (or  vice  versa).  It  is 
therefore  suggested  that  both  enhancements  be  implemented  into  the  RTNEPH 
coincidentally  so  as  to  minimize  the  length  and  complexity  of  the  readjustment 
periods  that  the  RTNEPH  will  undergo. 

At  first  glance  the  attenuation  correction  computations  outlined  in  this  study 
may  appear  too  complex  and  time-consuming  for  practical  RTNEPH  use.  However, 
these  models  are  fast,  efficient,  and  very  accurate.  The  most  intensive  transmittance 
computations  will  be  performed  on  a  whole  mesh  grid  basis.  For  example,  to 
compute  an  OLS-Thermal  band  transmittance  profile  and  subsequent  clear-column 
band  radiance  for  a  15-layer  atmosphere  takes  FASCODE  114  seconds  of  CPU  time 
on  a  3  MIPS  MicroVAX  3600  computer;  it  takes  the  multivariate  polynomial 
regression  transmittance  model  0.1  CPU  seconds,  faster  by  a  factor  of  1000. 

A  final  determination  will  be  made  on  the  accuracy  of  the  RTNEPH  cloud/no- 
cloud  decisions  using  the  current  and  the  proposed  attenuation  correction  procedures. 
Comparisons  of  cloud  fields  derived  using  each  procedure  will  be  made  directly 
against  interactively-derived  cloud  truth  fields  [see  Gustafson  and  Felde  (1988)  for  a 
description  of  how  RTNEPH  cloud  truth  fields  are  generated].  Using  these  fields, 
biases  inherent  in  the  proposed  attenuation  correction  procedures  can  be  discovered 
and  addressed.  It  is  believed  that  an  improved  attenuation  correction  model  will 
significantly  increase  the  quality  of  the  RTNEPH  IR  processor  clear/cloud  analysis. 
The  manually  intensive  process  of  updating  four  RTNEPH  attenuation  correction 
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and  bias  lookup  tables  will  be  reduced  to  at  most  one  bias  table,  with  the  other  three 
tables  being  replaced  explicitly  by  the  techniques  proposed  in  this  paper.  The 
proposed  attenuation  correction  techniques  can  be  easily  transitioned  to  new  sensors 
by  simply  changing  the  appropriate  sensor  response  function  values  fn  in  Eq.  (19c).  In 
the  same  way,  data  from  the  thermal  AVHRR  channels  at  3.7,  10.7,  and  11.8p,m  can 
be  incorporated  easily  into  the  RTNEPH  ER  processor  with  a  minimal  recoding 
effort. 
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